Berry-phase theory of polar discontinuities at oxide-oxide interfaces 
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In the framework of the modern theory of polarization, we rigorously establish the microscopic 
nature of the electric displacement field D. In particular, we show that the longitudinal component 
of D is preserved at a coherent and insulating interface. To motivate and elucidate our derivation, we 
use the example of LAO/STO interfaces and superlattices, where the validity of the above conserva- 
tion law is not immediately obvious. Our results generalize the "locality principle" of constrained-D 
density functional theory to the first-principles modeling of charge-mismatched systems. 
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In classical electrostatics, the normal component of the 
electric displacement field D is preserved, 



(D 1 -D 2 )-n = 0, 



(1) 



at an insulating (charge-carrier-free) planar interface be- 
tween two homogeneous insulators. Here Di and D 2 
are the values of the macroscopic electric displacement 
in material 1 and 2 and n is the plane orientation. 
While this is in principle a macroscopic law, recent 
first-principles calculations (e.g., on multicomponcnt per- 
ovskite superlattices )* have shown that it is applicable 
even at the microscopic level. Given this success, one 
would be tempted to adopt Eq. |T|) in full generality for 
the description of insulating interfaces at the nanoscale, 
where the conservation of D = D • n facilitates the inter- 
pretation and modeling of many phenomena. 

There are, however, a number of cases which have 
emerged recently where the applicability of Eq. |T]) is 
not immediately obvious. The prototypical example is 
that of a charge-mismatched interface between two crys- 
talline insulators. Polar interfaces have been the ob- 
ject of intense research in the past few years, motivated 
by the recent observation of two-dimensional metallic- 
ity at interfaces between LaA103 (LAO) and SrTiC>3 
(STO).-»2. Interestingly, first-principles calculations have 
recently demonstrated that the LAO/STO interface can 
remain insulating under certain conditions, and in such 
a regime the oxide lattice undergoes rather unusual re- 
laxation patterns. For example, in the case of a thin film 
of LAO deposited on the (001) surface of STO, strong 
"ferroelectric-like" polar distortions were reported in the 
LAO overlayer while the substrate remains essentially un- 
perturbed.— Conversely, in periodic LAO/STO superlat- 
tices it was shown that the LAO and STO components 
spontaneously polarize in opposite directions, the largest 
distortions occurring now in STO£ The source of this po- 
larization discontinuity AP is understood: it is induced 
by electric fields arising from extra interface charges ±e/2 
(for IV-III and II-III interfaces respectively) resulting 
from the "polar discontinuity" between the II-IV (STO) 
and III-III (LAO) perovskites. Similar issues arise for 
the case of II-IV/I-V (e.g., STO/KNb0 3 ) superlattices^ 



In these examples one is seemingly forced to conclude 
that Eq. (p} is violated, since the AP gives rise to a cor- 
responding discontinuity in D that is inconsistent with 
Eq. p. 

One way to resolve this issue is to accept a def- 
inition of "free charge" as including everything other 
than bound charge, as suggested by some authors (e.g., 
Ref. 0, Sec. 4.3.1). Then Eq. © is fixed by adding a 
"free charge" term on the right to represent the polar- 
discontinuity charge, even though this is anything but 
"free." Another possibility, proposed by Murray and 
Vanderbilt pi is to introduce a compositional charge den- 
sity Pcomp, which is distinct both from bound charge 
(Pbound = — V • P) and from truly free charge (associ- 
ated with charge carriers), and write 



Ptot — Pfrcc "T /Abound t /?comp- 



(2) 



With V • D = 47r(pf rec + Pcomp), Eq. ^ becomes, for an 
insulating interface, 



(Di - D a ) • n = Ana 



comp ■ 



(3) 



While these approaches already provide a satisfactory 
account of the phenomena described in Refs. HUB there 
are good reasons to seek an alternative description. First, 
both of the above approaches are somewhat awkward, ei- 
ther introducing a third kind of charge, or describing as 
"free" a form of charge that is fixed. Second, it is nec- 
essary to assess the range of applicability of Eq. (Q]), i.e. 
to identify a general criterion for deciding in what cases 
it should be replaced by Eq. ([3]). Third, one would like 
to establish a true microscopic definition of all quantities 
appearing in Eq. @, particularly p comp (r) and Abound (r). 
Given the ever-increasing interest in perovskite thin films 
and superlattices, such an analysis is particularly urgent 
in order to provide timely support for the experimental 
work with appropriate modeling tools. 

Here we show that the incompatibility between Eq. (Q]) 
and the "polar discontinuity" can instead be elegantly re- 
solved in the framework of the Berry-phase modern the- 
ory of polarization. In particular, we show that Eq. ([T]) 



2 



is exact without the need for extensions, once D is ex- 
pressed in terms of the formal (rather than the effective) 
macroscopic polarization, and is a direct consequence 
of the interface theorem^ This result puts the "locality 
principle" of constrained-D density functional theory^ on 
a firm theoretical basis, and generalizes its scope to the 
case of charge-mismatched superlattices and heterostruc- 
tures. In the following, we will first introduce some basic 
properties of the formal polarization in periodic insula- 
tors. Based on these properties, we will then establish the 
link between Eq. JT]) and the interface theorem. Finally, 
we will demonstrate these ideas in practice by performing 
calculations of a model LAO/STO superlattice. 

In an independent-electron (or mean-field) picture, the 
total charge density of an insulator can be expressed as 
a sum of contributions from ion cores and Wannier func- 
tions. Let r a and represent the ion core and Wannier 
center locations, respectively, for a choice of "basis" that 
tiles the crystal under primitive translations a l5 a 2 and 
a3. The dipole moment of this discrete set of charged 
objects divided by the cell volume defines the macro- 
scopic polarization Pi&Iii We shall be concerned with a 
single component of the polarization, i.e., its projection 
in the direction n parallel to a 2 x a3; let this be direction 
x. Then we have an essentially one-dimensional problem 
with P = P x given by 
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(4) 



where Q a and — e are the charges of a given ionic core and 
Wannier function, respectively, and Q is the cell volume. 
While this is rigorous, it has a degree of arbitrariness 
in that one may choose to include in the crystal basis 
any of the infinitely repeated images of each Wannier 
function or ion core. As shown schematically in Fig.QJa), 
this implies that P is a multi-valued function; it is only 
defined modulo a "quantum of polarization" AP = e/S, 
where S = |a 2 x ag| is the cell surface area. 

The polarization in Eq. (j4|) is the formal polarization 
corresponding to the raw result of a Berry-phase calcu- 
lation (in contrast to the effective polarization, defined 
relative to a centrosymmetric reference) ^ It is essential 
to understand that the formal polarization does not nec- 
essarily vanish in a centrosymmetric material (i.e., P=0 
is not contained in the lattice of allowed values). Practi- 
cal realizations of this situation are III-III perovskites like 
LAO in their cubic five-atom reference structure. With 
x along the (100) direction, the lattice of allowed values 
is P = ±e/2S, ±3e/2S, as shown schematically in 
Fig. [T](b) . This occurs essentially because the individual 
LaO and A10 2 layers have formal charges of ±e. In con- 
trast, II-IV perovskites like STO, which have formally 
neutral AO and B0 2 layers, have P = 0, ±e/S, as 
illustrated in Fig. [IJa). We now show that this funda- 
mental difference in the respective P-lattices of LAO and 
STO is the key to reconciling Eq. (fTJ) and the "misbehav- 
ior" of charge-mismatched interfaces. 

Consider a coherent insulating 12 interface between 
LAO and STO with perfect (1 x 1) periodicity, i.e. a 2 and 
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FIG. 1: (Color online) Schematic representation of formal 
polarization in the bulk (a-b) and at interfaces (c-d). Blue 
and red vertical lines represent BO2 and AO layers in a II-IV 
perovskite (a); green and brown are the same but in a III- 
III perovskite (b). Black circles represent electronic Wannier 
centers (not all are explicitly shown). Different choices of 
basis (grouping) lead to values of P differing by AP = e/S. 
(c) Choice of basis such that no charges are left in the interface 
region, (d) Choice of bias leaving net-neutral interface region 
(grey). 



a3 are common to both crystals and lie on the interface 
plane. Fig.[TJc-d) shows an example with a Ti0 2 :LaO in- 
terface termination, but similar considerations apply for 
other terminations. Then it is always possible to group 
ions and Wannier functions into a basis for STO and one 
for LAO (yellow areas in the figure) such that the left- 
over interface region (grey shaded area) is overall charge 
neutral; two examples are provided in the figure. We 
use these basis conventions to define the bulk polariza- 
tion in STO (Pi) and LAO (P 2 ) via Eq. gj). Because of 
the neutrality of the interface region, it follows that the 
macroscopic density of bound charge at the interface is 



0~ba 



Pi -P 2 - 



(5) 



This is essentially the "interface theorem" of Ref. [8|, a 
central result of the modern theory of polarization. 

Eq. ([S]) implies that the bound charge at the interface 
is determined by the discontinuity of formal bulk po- 
larizations of the two participating materials, removing 
the need to complement the theory with p comp £L At the 
Ti0 2 :LaO interface sketched in Fig. QJ c-d), for example, 
Eq. ([5]) would yield a abound of exactly half an electron 
per unit cell if both bulks were centrosymmetric, con- 
sistent with the heuristic arguments of the "polar catas- 
trophe" 3 model. However, macroscopic electric fields £\ 
and £2 will typically be present in materials 1 and 2, since 
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Gauss's theorem implies that 



-£i + £2 



(6) 



and the self-consistent solution is the one in which the 
polarizations Pi and P2 are the equilibrium values in the 
corresponding fields £\ and £2- These polar distortions, 
in turn, screen the discontinuity in P. Introducing the 
electric displacement^ 



D = £ + 4ttP , 



(7) 



it follows from Eq. ([6|) that D\ — D2 = 0, which is ex- 
actly Eq. (JXJ) . [Note that the multivalued character of P 
propagates to the D through Eq. ([7]) ; D is lattice- valued 
with a "quantum" AD — e/AirS.] Thus, the phenomena 
discussed in the introduction emerge as a consequence to 
Eq. Jl}, rather than as an exception to it. 

To substantiate this interpretation, where all charge 
not associated with free carriers is counted as bound 
charge, pbound(r) = ptot(r), it is useful to look back at 
the theory developed in Ref. [l5|. The authors defined a 
"local" polarization through 16 



dP(x) 
dx 



Pbound (x*) , 



(8) 



where the bar indicates planar averaging over the yz 
planes and "macroscopic" averaging^ along the het- 
erostructure stacking direction x. This leads to a mi- 
croscopic definition of the displacement field, D(x) = 
£(x) + 4nP(x), where the local electric field is given by 
the microscopic Maxwell equations, 



d£{x) 
dx 



-^Ptot(x). 



(9) 



Combining Eqs. (JSJ) and © we obtain that D(x) must 
be constant throughout the superlattice, consistent with 
Eq. {l]). It is a l so straightforward to show that P(x) 
converges to P\ or P2 sufficiently far from the interface. 

These arguments, therefore, generalize both the "lo- 
cality principle" of constrained-D theory^ and the the- 
ory of local dielectric response^ to the case of charge- 
mismatched systems. However, this result comes at a 
price: One must accept that the formal polarization in 
a centrosymmetric perovskite like LAO does not vanish 
(and might even be non-zero in centrosymmetric STO, 
depending on the branch choice). While this may be 
counterintuitive, it is an established aspect of the mod- 
ern theory of polarization^ and has been already crucial 
for answering important questions in the physics of com- 
plex oxides such as BiFeOsJ^ Here we show that the 
"half-quantum" nature of a III-III (or I-V) perovskite is 
not merely a technical annoyance; rather, it acquires a 
well-defined physical meaning through Eq. ([5]). 

To elucidate in practice the continuity of the formal 
polarization at a polar interface, we now consider explic- 
itly a periodic superlattice composed of four layers each of 
STO and LAO in a tetragonal lxl supercell with the in- 
plane lattice parameter set to the theoretical equilibrium 
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FIG. 2: (Color online) Top panel: cation-oxygen rumplings, 
Sao and <5bo2, in a given layer of the (LAO)4/(STO)4 su- 
perlattice. The red boxes indicate the grouping of the layers 
adopted for defining Pj and Pj. Bottom panel: cell polariza- 
tions Pj (open symbols), Pj (filled symbols) and macroscopic 
value PBerry (dashed horizontal line). 



value of ao = 7.274 a.u. for cubic STO. Our calculations 
are performed within the local-density approximation of 
density-functional theory and the projector-augmented- 
wave method^ with a plane-wave cutoff of 60 Ry. The 
Brillouin zone is sampled with a (6 x 6 x 1) Monkhorst- 
Pack grid. We relax the ionic forces and the out-of-plane 
stress to a tolerance of 10 -5 Ha/bohr and 10 -7 Ha/bohr 2 , 
respectively. After relaxation the system remains insulat- 
ing with a large gap of ~ 1.5 eV, in qualitative agreement 
with Ref. 0. We then inspect the local polarization profile 
using two contrasting approaches as follows. 

First, we define a local formal polarization Pj for each 
of the eight cells j by inserting the coordinates of the 
ionic positions and Wannier centers of the cell in question 
into Eq. Qi 2 ^ The definition of a "cell" is fixed by the 
previous choice of basis for each bulk material. The cal- 
culated Pj values are shown as the filled circles in Fig. [21 
where it is assumed that the chosen basis was the one 
sketched in the top panel of Fig. [21 also corresponding 
to Fig. [ITc). The total Berry-phase polarization of the 
supercell, PBerry = 1/8 ^2j Pj, is then simply the average 
of the local Pj values. (Our definition of Pj is closely re- 
lated to the prescription of Ref. [2l|, except that a further 
subdivision into contributions from AO and BO2 layers 
was introduced there.) 

Second, by using the same grouping convention, we 
define the local effective polarization Pj of cell j as 



P, 



ZlAx c 



(10) 



where Z* and Ax a are the Born effective charge and 
atomic displacement of atom a respectively. (Displace- 
ments are defined relative to a stack of ideal centrosym- 
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metric cells.) The resulting Pj values are plotted as the 
open symbols in Fig O 

Strikingly, the Wannier-based formal polarization Pj 
is rather uniform, with little departure from the average 
value Pecrry = —0.365 C/m 2 shown as the dashed line 
in Fig. [2] This confirms the validity of our arguments: 
At electrostatic equilibrium the system tries to minimize 
the macroscopic electric fields present on either side of 
the interface by making the discontinuity in the "formal" 
P as small as possible. Interestingly, in our heterostruc- 
ture, the LAO and STO components acquire oppositely 
oriented structural distortions (upper panel of Fig. [5]) in- 
order to achieve this goal. This is reflected by the abrupt 
discontinuity in the "effective" local polarization pattern 
Pj . Note that in the STO region Pj and Pj almost coin- 
cide, since our branch choice for the formal polarizations 
implies P = P = for bulk STO in its centrosymmet- 
ric ground state. The difference between Pj and Pj in 
the LAO region of the plot (Fig. \2§ corresponds to half a 
quantum of polarization, as expected. 

To corroborate our conclusions, we now show that the 
local properties deep in the LAO and STO regions are 
determined by the macroscopic D, and do not directly 
depend on the specific details of the interface. To that 
end, we calculate (within the same computational param- 
eters and symmetry constraints) the equilibrium struc- 
ture of bulk LAO and bulk STO with D constrained 14 
to take the value D=— 0.365 C/m 2 extracted from the 
heterostructure in electrostatic equilibrium. The av- 
erage bucklings 5LaO=0-078A and <5aio2 = 0.044A ex- 
tracted from cells 6-7 are in excellent agreement with 



bulk LAO values of 0.078 A and 0.043 A respectively, 
while (5 Sr o=-0.140 A and <5 TiO 2=-0.140 A (coinciden- 
tally the same for both layers) extracted from cells 2-3 
similarly match the bulk STO values of —0.137 A and 
—0.140 A. The largest discrepancy is ~3mA, confirming 
the soundness of our arguments. A similar reasoning 
can be used to interpret the results of Refs. and [22c 
when D=0 is enforced by symmetry, an electric field 2 ^ of 
£=0.24 V/A and significant ferroelectric-like rumplings 4 
(<5LaO=0.26 A, <5aio2=0.15 A) were reported in the LAO 
layer. To check whether these results are consistent 
with our arguments, we repeated our bulk LAO calcu- 
lations at D=0 and found £=0.238 V/A, <5 LaO =0.24lA 
and <5aiO2=0.143 A, in excellent agreement with the lit- 
erature values. 

In conclusion, we have demonstrated that the conser- 
vation of the longitudinal component of D expressed in 
Eq. ||T|) is a microscopic law that explains and predicts 
the behavior of insulating interfaces and superlattices in 
a variety of electrical boundary conditions. This fun- 
damental principle, together with state-of-the art finite 
field approaches is a powerful theoretical tool to com- 
plement and guide experiments in the emerging field of 
interface nanoengineering. Our ideas are very general, 
and can readily be applied to a wide range of physical 
systems such as, e.g., nitride superlattices where the is- 
sue of interface polarity poses important technological 
challenges^ 
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